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Abstract. The study of shell-type supernova rem- 
nants is a key science focus for the VERITAS TeV 
telescope array. Supernova remnants (SNRs) are 
widely considered to be the strongest candidate for 
the source of cosmic rays below the knee around 
10 15 eV. This presentation will highlight new VER- 
ITAS results including new measurements of the 
spectra of Cas A and IC 443. These results and 
their implications for the nature of the cosmic rays - 
hadronic or electronic - accelerated in the remnants 
will be discussed. 
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I. Introduction 

Fermi acceleration in supernova remnants (SNRs) has 
long been considered as one of the likely explanations 
for the origin of cosmic rays (CRs) up to the knee at 3 • 
10 15 eV. This scenario drives the observations of SNRs 
in the very-high-energy (VHE, E > 100 GeV) gamma- 
ray band. Measurements of the VHE morphology and 
energy spectrum, when combined with multiwavelength 
observations, hold promise to constrain models of diffu- 
sive shock acceleration and the potential to discriminate 
between gamma-ray production by leptonic and hadronic 
cosmic rays. For example, the deep H.E.S.S. observa- 
tions of RX J17 13.7-3946 reveal an energy spectrum 
extending beyond 30 TeV, implying particle acceleration 
to well beyond 100 TeV Q. Unfortunately, in many 
cases (including RX J17 13.7-3946) uncertainty in the 
environment (age and distance of remnant, density of 
interstellar medium) makes it difficult to discriminate 
conclusively between hadronic and leptonic scenarios for 
the gamma-ray emission. 

This paper discusses detections of VHE gamma rays 
from the supernova remnants Cassiopeia A and IC 443. 
An upper limit on the "Forbidden Velocity Wing" FVW 
190.2+1.1, a possible supernova remnant, is discussed 
in |27|. A summary of the full VERITAS observation 
program on supernova remnants and pulsar wind nebulae 
is given in OTIl in these proceedings. 

II. Observations and Analysis 

VERITAS 1 26] consists of four 12-m telescopes lo- 
cated at an altitude of 1268 m a.s.l. at the Fred Lawrence 
Whipple Observatory in southern Arizona, USA (31° 40' 
30" N, 110° 57' 07" W). Each telescope is equipped 
with a 499-pixel camera of 3.5° field of view. The array, 
completed in the spring of 2007, is sensitive to a point 



source of 1% of the steady Crab Nebula flux above 
300 GeV at ha in less than 50 hours at 20° zenith angle. 

Observations are taken in wobble mode fl9|, with the 
telescope pointing offset by 0.5°-0.7° from the target 
direction for roughly equal amounts of time in each of 
four equally spaced directions. 

The data analysis is described more fully in ifTTl . 
In brief, pixel gains are flatfielded using laser flashes 
recorded during nightly calibration runs. A two-pass 
cleaning algorithm with thresholds of 5 and 2.5 times the 
pedestal RMS is used; isolated pixels which survive the 
image cleaning are then eliminated. The images are pa- 
rameterized according to their first and second moments 
11251 . and then the shower direction and core location 
are determined. Lookup tables provide expected values 
for the image width, length, and energy, and are used 
to construct the gamma-hadron separation parameters 
Mean Scaled Width (MSW) and Length (MSL), and an 
estimate of the shower energy [ 33 1 . Images are included 
in the reconstruction if they have at least five pixels, if 
the centroid is less than 1.43° from the center of the 
camera, and if the image size (sum of the pixel charges) 
is at least 75 p.e. At least two good images per event 
are required. Events in which only the two telescopes 
sharing the shortest baseline have good images are 
rejected, as their 35-m separation is too small for good 
stereo reconstruction. Cuts on MSW and MSL have been 
optimized for weak sources with a Crab-like spectrum. 
The background is estimated using the "ring background 
model" [5|. Significances are calculated according to 
equation 17 of l35l . All results are verified with an 
independent analysis. 

III. Cassiopeia A 

Cassiopeia A (Cas A) is a young remnant, believed 
to be about 300 years old and at a distance of 3.4 kpc. 
It is the brightest and among the best-studied radio 
sources in the sky (32). The synchrotron radiation from 
Cas A extends from radio wavelengths to hard X-rays 
IfTTl . |T8). It has a diameter of 5 arcmin, making it a 
point source for TeV telescopes. Cas A was not detected 
by EGRET, and is not seen in the Fermi Gamma- 
ray Space Telescope (FGST) bright-source list (TJ. Cas 
A was first detected in TeV gamma rays by HEGRA 
|4|, with a flux of ~ 3.3% Crab above 1 TeV. The 
detection was confirmed by MAGIC |9| at 5.3 a in 
a 47-hour observation, yielding a flux above 1 TeV 
of (7.3 ± 0.7 sta t ± 2.2 sas ) x 10" 13 cm^s" 1 and a 
spectral index of 2.3 ± 0.2 sta t ± 0.2 sys , consistent with 
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Fig. 1. Cas A significance map, showing a strong point-like signal 
at the location of Cas A. 
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Fig. 2. Cas A 8 2 distribution, indicating the number of events as a 
function of the square of the separation angle from the center of Cas 
A. The filled circles are for the on-source region and the empty circles 
are an estimate of the background via the reflected-region method. 



the HEGRA results. Cas A is an interesting TeV SNR 
because of its low age, detailed observations across all 
wavelengths, and relatively simple local environment, 
making it attractive for constraining models of particle 
acceleration in shocks. 

Cas A was observed in October/November, 2007 with 
four telescopes. After data quality selection (for rate 
stability and good weather) and correction for deadtime, 
an exposure of 20.3 hours remains with an average 
zenith angle of 30°. 

Preliminary results, previously show in [30], are given 
here, and final results will be presented at the ICRC. Fig- 
ure Q] shows the significance map for the Cas A field and 
Figure [2] shows the distribution of events as a function 
of the square of the angular distance from the source 
location. Clear peaks are seen at the location of Cas 
A, with a significance of 9.8 a and a flux (determined 
by comparing to the rate of gamma rays from the Crab 
Nebula) of ~ 3 % Crab, consistent with the HEGRA 
and MAGIC results. The emission is consistent with a 
point source within the instrument's angular resolution, 
and the spectral analysis is in progress. 

IV. IC443 

IC 443 is one of the classic examples of a SNR 
interacting with a molecular cloud, and this interaction 
greatly affects its evolution. VHE emission coincident 
with the site of the shock/cloud interaction arguably 
provides a smoking gun for the acceleration of hadronic 
cosmic rays. IC 443 is at a distance of 1.5 kpc and 
has a diameter of ~ 0.75°. The northeast hemisphere 



is interacting with an Hll region, and the southwest 
hemisphere is interacting with a molecular cloud at 
several points. The cloud has a mass of ~ 10 4 M Q 
[Q3], mostly along the line of site to IC 443; maser 
emission has been observed in this region lfl6l . The age 
of IC 443 is unclear; various estimates place it in the 
range 3-30 kyr (see, e.g., ED, EQ). The PWN CXOU 
J06 1705. 3+222 127 is located at the southern edge of the 
remnant J37), 03, ED- Both MAGIC and VERITAS 
have reported TeV emission from the central region of 
the remnant, coincident with the maser emission J8), 
). An EGRET GeV source also overlaps the remnant 
and recent detections by AGILE (59) and FGST 
ID have been reported. 

IC 443 was observed by VERITAS during two pe- 
riods: with three telescopes in February /March, 2007 
(during the commissioning phase of the array) for 20 
hours; the PWN was the target for these observations. 
A further 25 hours of observations were taken in Octo- 
ber/November, 2007 using the complete four-telescope 
array and targeting the location of the emission observed 
previously, 06 16.9 +22 33. After quality selection and 
deadtime correction, the data set has a livetime of 37.1 
hours and an average zenith angle of 18°. 

One systematic issue particular to this data set is 
the presence of two bright stars in the IC 443 field: 
Eta Gem (V band magnitude 3.31), located 0.53° from 
the PWN; and Mu Gem (V band magnitude 2.87), 
located 1.37° from the PWN. The high flux of optical 
photons from the stars requires that several PMTs in 
each camera be turned off during the observations and 
produces increased noise levels in the signals of nearby 
PMTs. These effects reduce the exposure in the vicinity 
of the stars and degrade the angular resolution. To 
counteract these effects while studying the morphology 
of the source, a high telescope-multiplicity requirement, 
described below, was selected. 

A standard point-source analysis of the full data set 
yields a maximum significance of 8.3 standard devia- 
tions (a) before trials (7.5 a post trials, accounting for 
a blind search over the region enclosed by the shell of 
IC 443) at 06 h 16 m 49 s +22°28'30" (J2000) and a signif- 
icance of 6.8 a at the location of MAGIC J06 16+225 
(see (J) for full details). Table J] summarizes the results 
of this analysis at the location of maximum significance, 
listing the counts falling within the point-source inte- 
gration radius of 0.112° (on), the counts integrated in a 
background ring spanning radii 0.6-0.8° (off), the ratio 
of on exposure to off exposure (a), and the resulting 
number of excess counts and significance. Figure [3] 
shows the inner 0.8° of the excess map. To study 
the source morphology, the best-reconstructed gamma 
rays are selected by requiring that an event must have 
images from all participating telescopes survive the pre- 
selection cuts. This requirement increases the analysis 
energy threshold by ~15% and reduces the excess by 
~40%. The centroid and intrinsic extension of the excess 
are characterized by fitting an azimuthally symmetric 
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two-dimensional Gaussian, convolved with the PSF0 of 
the instrument, to an acceptance-corrected uncorrelated 
excess map with a bin size of 0.05°. The PSF has a 68% 
containment radius of 0.11°. The centroid determined in 
this way is located at 06 h 16 m 51 s +22°30'll" (J2000) ± 
0.03° sta t ± 0.08° sys , consistent with the MAGIC posi- 
tion, and the extension is 0.16° ± 0.03° stat ± 0.04° sys . 
Due to the effects of the nearby bright star and the non- 
Gaussian nature of the excess, a large systematic error 
is currently assigned to the centroid. The difference in 
extension between this work and the point-like source 
detected by |8) can be explained by the difference in 
angular resolution and sensitivity between VERITAS 
and MAGIC, the latter of which has an angular reso- 
lution of ^0.15° for 68% containment and a sensitivity 
above 200 GeV of 2% of the Crab Nebula flux in 50 
hours [ 1 1 . The photon spectrum, integrated within a 
radius of 0.235°, is shown in Figure [4] The photon 
spectrum is well fit (x 2 /ndf = 3.1/3) by a power law 
dN/dE = N x (£'/TeV)" r with a normalization of 
(8.38 ± 2.10 stat ± 2.50 sys ) x 10~ 13 TeV" 1 cm" 2 s" 1 
and an index of 2.99 ± 0.38 s t a t ± 0.30 sas . The in- 
tegral flux above a threshold energy of 300 GeV is 
(4.63 ± 0.90 stat ± 0.93 sys ) x 10~ 12 cm" 2 s" 1 (3.2% 
of the Crab Nebula flux). 

TABLE I 

Analysis results for VERITAS observations of IC 443. 





on 


off 


a 


excess 


significance 

(<y~) 


data set 1 


393 


5152 


0.0610 


78.7 


4.1 


data set 2 


609 


7331 


0.0601 


168.1 


7.3 


combined 


1002 


12483 


0.0604 


247.5 


8.3 



The TeV centroid is ~ 0.15° from the position of 
the PWN and ~ 0.03° from the nearby, bright maser 
emission coincident with clump G of ESI . Weaker 
regions of maser emission lie along the southern rim 
of the remnant, east of the PWN |24|. The VHE emis- 
sion coincides with the densest region of the molecular 
cloud. While the Fermi source 0FGL J06 17. 4+2234 is 
displaced from the centroid of the VHE emission by 
^0.15°, this is consistent at the 95% level with the 
combined errors quoted by Fermi and VERITAS. 

The VHE gamma-ray emission observed by VER- 
ITAS and MAGIC is offset from the location of the 
PWN by ~ 10-20 arcmin, similar to other offset PWNe 
such as HESS J 1825- 137 |6). The emission is consistent 
with a scenario in which the VHE emission arises from 
inverse Compton scattering off electrons accelerated 
early in the PWN's life, as in |Q2|. If the VHE emis- 
sion is indeed associated with IC 443 at a distance of 
1.5 kpc, the luminosity in the energy band 0.3-2.0 TeV 
is 4 x 10 32 erg s -1 . The spin-down luminosity of the 
pulsar has been estimated as ~10 36 erg s _1 11371 . |[T4l 

'The PSF is characterized as a sum of two, two-dimensional 
Gaussians describing a narrow core and a broader tail. It is determined 
from data taken on the Crab Nebula, which is a point source at these 
energies |10|. 
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Fig. 3. Inner 0.8° of the acceptance-corrected excess map for the 
IC 443 field. The black cross-hair indicates the centroid position and 
its uncertainty (statistical and systematic added in quadrature), and 
the white cross-hair likewise indicates the position and uncertainty 
of MAGIC J0616+225 (8). Red contours: optical intensity (36). 
Thick black contours: CO survey |'28J; black star: PWN CXOU 
J061705. 3+222127 [37|; open blue circle: 95% confidence radius of 
0FGL J0617.4+2234; and filled black triangles: locations of OH maser 
emission f 1161 . 1241 . J. W. Hewitt, private communication). The white 
circle indicates the PSF of the VERITAS array. 



and 5 x 10 37 erg s -1 (20|. If the PWN association is 
correct, the VHE luminosity in the 0.3-2.0 TeV band is 
less than ^0.04% of the spin-down luminosity, at the 
lower end of the ^0.01-10% range of observed VHE 
efficiencies for other PWNe [22|. 

Alternatively, Figure [3] can be interpreted as hadronic 
cosmic-ray acceleration and subsequent interaction with 
the molecular cloud, which would provide a high density 
of target material for the production of VHE gamma 
rays. The steep VHE spectrum can be explained either 
as an energy-dependent rate of diffusion of cosmic rays 
out of the cloud [3 j, or a low maximum energy to which 
particles were accelerated prior to the shock hitting the 
cloud. 

Further GeV and TeV gamma-ray observations will 
aid in elucidating the nature of the particle acceler- 
ation associated with IC 443. The ^0.15° separation 
between the centroids of 0FGL J06 17. 4+2234 and VER 
J0616. 9+2230, though not statistically significant, may 
hint at an energy-dependent morphology. This will be 
explored with future observations by VERITAS, and 
has the potential to discriminate between the leptonic 
and hadronic scenarios sketched above. The relocation 
of one of the telescopes during the summer of 2009 
as part of the VERITAS upgrade program 11381 will 
improve both the angular resolution and sensitivity of 
the array, assisting with these studies. Whether or not 
0FGL J06 17. 4+2234 and VER J06 16.9+2230 are di- 
rectly related, a break in the energy spectrum between 
the GeV and TeV bands is required. Future observa- 
tions by VERITAS will improve the statistics above 
100 GeV, anchoring the measurement of the VHE cut- 
off to the gamma-ray emission. The combined spectrum 
from 100 MeV to several TeV will provide powerful 
constraints on the emission mechanism(s) [21 j. 
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Fig. 4. Spectrum of IC 443, scaled by i? 2 . Red points are the 
VERITAS spectrum (red tick marks along top indicate bin edges). The 
red line is a power-law fit (see text), with residuals to the fit plotted 
in the lower box. The MAGIC spectrum is indicated by the gray band 
and extends down to 90 GeV. VERITAS error bars and MAGIC error 
band reflect statistical errors only. The Crab spectrum (V. Acciari et 
al. 2009, in preparation) is shown as a dashed line for comparison. 

V. Conclusions 

VERITAS has made strong detections of both the 
Cassiopeia A and IC 443 supernova remnants. Cas A 
is observed to have a flux of ~ 3% Crab and appears 
pointlike within the instrument's resolution. These re- 
sults are consistent with previous observations, and final 
results, including a spectrum and an upper limit on the 
source extension, are in preparation. 

The emission from IC 443 is extended, yielding a 
sigma of ~ 0.16° when fit with a two-dimensional 
gaussian. It has an integral flux of 3.2% Crab above 
300 GeV. The TeV emission is coincident with both 
the densest part of the molecular cloud that IC 443 is 
interacting with and the maser emission observed within 
that cloud, strongly suggesting that the gamma rays are 
produced by the interaction of hadronic cosmic rays with 
cloud material. An alternative explanation in terms of an 
association with the PWN CXOU J061705.3+222127 is 
also possible. 

The VERITAS survey of the Cygnus region of the 
Galactic plane, spanning Galactic longitude 67 < I < 82 
and latitude — 1 < b < 4, covers a number of other 
supernova remnants such as 7Cygni and CTB 87 and is 
discussed elsewhere in these proceedings |42|. 
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